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ABSTRACT 

The second moment of the H/3 emission line is calculated for 329 narrow line Seyfert 
1 galaxies (NLSls) selected from the Sloan Digital Sky Survey (SDSS), which is used 
to calculate the central supermassive black hole (SMBHs) mass of each. We find that 
the second moment depends strongly on the broader component of the H/3 line profile. 
We find that for the NLSls requiring two Gaussians to fit the H/3 line the mean value 
of the SMBH mass from the H/3 second moment is larger by about 0.50 dex than that 
from the full width at half maximum (FWHM) . Using the gas velocity dispersion of the 
core/narrow component of [O III] A 5007 to estimate the stellar velocity dispersion, cr», 
the new mass makes NLSls fall very close to the Mbh — 0* relation for normal AGNs. 
By using er* measured directly from SDSS spectra with a simple stellar population 
synthesis method, we find that for NLSls with mass lower than 10 7 M Q , they fall only 
marginally below the Mbh — o* relation considering the large scatter in the mass 
calculation. 
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1 INTRODUCTION 

Narrow line Seyfert 1 galaxies (NLSls) are thought to be 
a special subclass of active galactic nuclei (AGNs) harbor- 
ing relatively small but growing supermassive black holes 
(Mbh , SMBHs), compared to other broad line Seyfert 1 
galaxies (BLSls; e.g., Osterbrock & Pogge 1985; Boiler et 
al. 1996; Mathur et al. 2000). Whether NLSls follow the 
well-known Mbh — cr* (or Mbh — Ltmige) relation defined in 
inactive galaxies is a question open to debate, where the stel- 
lar velocity dispersion (<t„) is measured at an eighth of the 
effective radius of the galaxies (e.g., Tremaine et al. 2002 
Mathur Kuraszkiewicz & Czerney 2001; Bian & Zhao 2004 
Grupe & Mathur 2004; Botte et al. 2005; Barth et al. 2005 
Watson Mathur & Grupe 2007; Ryan et al. 2007; Komossa 
& Xu 2007). It remains also a question for other types of 
AGNs (e.g., Nelson et al. 2000; Greene & Ho 2006; Shen et 
al. 2008; Woo et al. 2008). In investigating the Mbh - cr» 
relation for NLSls or/and other AGNs, the method of de- 
terminging Mbh and a* is very important. 

The central SMBH mass in AGNs is a key parameter to 
understand the nuclear energy mechanism as well as the cos- 
mic formation and evolution of SMBHs and their host galax- 
ies (e.g. Rees 1984; Gebhardt et al. 2000; Ferrarese & Merrit 
2000; Tremaine et al. 2002). In the past two decades, there 



has been striking progress in finding more reliable methods 
to calculate SMBHs masses in AGN through the line width, 
AV, of H/3 (or Ha, Mg n , C IV ) from the broad line re- 
gion (BLR) and the BLR size, -Rblr (e.g., Kaspi et al. 2000; 
McLure & Dunlop 2004; Bian & Zhao 2004; Peterson et al. 
2004; Greene & Ho 2005b). 

Introducing the scaling factor, /, to characterize the un- 
clear kinematics and the geometry of the BLRs, the SMBH 
mass is calculated by: 



Mbh = / 



AblrAU 2 
G 



(1) 



The uncertainties of SMBHs masses in AGN from Equation 
1 are mainly from the uncertainties in Rblr, f and AV. 

Much effort has been focused on determining -Rblr, 
from the reverberation mapping method or empirical size- 
luminosity relations (e.g., Kaspi et al. 2000; 2005; Vester- 
gaard & Peterson 2006; Bentz et al. 2006). There are mainly 
two ways to parameterize the line widths of broad emis- 
sion lines, i.e., the full width at half maximum (FWHM) 
and the second moment (ai ine ). For a Gaussian line profile, 
FWHM/(7ii Ilc = V81n2 « 2.35; while for a Lorentzian pro- 
file, aiine — ► oo. However, people usually find that one Gaus- 
sian component provides a poor fit to the H/3 emission line 
profile after subtracting the contribution from narrow line 
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region (NLR), especially for NLSls (e.g., Rodriguez- Ardila 
et al. 2000; Dietrich et al. 2005; Mullaney & Ward 2008). 
Salviander et al. (2007) used a Gauss-Hermite function to 
measure the H/3 FWHM (also see McGill et al. 2008). Net- 
zer & Trakhtenbrot (2007) measure FWHM from the two- 
Gaussian fits to the H/3 profile (see also Mullaney & Ward 
2008). For 12 NLSls, Dietrich et al. (2005) suggested that 
their broad emission line profiles are well represented by 
employing two Gaussian components, while the Lorentzian 
profile does not do as well for the core and wing of the H/3 
profile simultaneously. Based on the analysis of reverbera- 
tion mapping data, it was suggested that an ne rather than 
FWHM be used to characterize the line width (Fromerth & 
Melia 2000; Krolik et al. 2001; Peterson et al. 2004). 

For 16 AGNs with BLRs sizes from the reverberation 
mapping and possessing reliable a» measurements, Onken 
et al. (2004) determined the scaling factor, /, to make the 
reverberation-based SMBH masses consistent with the well- 
known Mbh — c* relation of inactive galaxies (Tremaine et 
al. 2002). They found that / is 5.5 ± 1.8 when <7ii ne from 
root-mean-square (rms) spectra is adopted as AV . Collin 
et al. (2006) proposed different / factors for emission lines 
with different ratios of FWHM to crime- They suggested 
/ = 3.85 ± 1.15 when o"ij llc from mean spectra is adopted as 
AV. It is also often assumed that the BLR gas has random 
orbits. Netzer (1990) suggested that / = 3 when FWHM/2 
is adopted as AV (e.g., Kaspi et al. 2000; 2005; Greene & 
Ho 2006). 

Although <r» is difficult to measure for AGNs because 
the nuclei outshine their hosts, there are a larger number 
of AGNs from the Sloan Digital Sky Survey (SDSS) with 
obvious stellar absorption features within 3" aperture spec- 
tra, which can be used to measure a* (e.g., Green & Ho 
2006, Shen et al. 2008). About the method to measure 
please refer to Bian et al. (2007) and the reference therein. 
The gaseous velocity dispersion (e.g., ofoTn]! °"[n n]) is of- 
ten used to as a proxy for <r* (e.g., Nelson & whittle 1996; 
Greene & Ho 2005a). 

In our previous work based on NLSls selected from the 
Sloan Digital Sky Survey early data release (SDSS EDR), us- 
ing H/3 FWHM to calculate the mass and [O ill] narrow/core 
FWHM to give the cr* , we found that the SMBH masses of 
NLSls deviated significantly from the well-known Mbh — cr„ 
relation (Bian & Zhao 2004, Bian et al. 2006). Here, we 
use the second moment of the broad H/3 profile to re- 
investigate SMBHs masses in NLSls. All of the cosmological 
calculations in this paper assume Ho = 70 km s~ Mpc~ , 
M = 0.3, and Ov = 0.7. 



2 SAMPLE AND DATA ANALYSIS 

We use a largest sample of about 2000 NLSls from SDSS 
DR3 (Zhou et al. 2006). Zhou et al. (2006) presented this 
sample of about 2000 NLSls selected from objects assigned 
as "QSOs" and "galaxies" in the spectroscopic database of 
SDSS DR3. The only criterion is that the broad component 
of H/3 or Ha is detected and is narrower than 2200 km s _1 
in FWHM. Using a Lorentzian function to model the broad 
H/3 profile, Zhou et al. (2006) obtained the FWHM and cal- 
culated SMBHs masses from the J?blr — £5100 relation of 



Kaspi et al. (2000). We use the latest SDSS spectra from 
data release 6 (DR6). 

We briefly outline our steps to do the SDSS spec- 
tral analysis. (1) We use simple stellar population (SSP) 
synthesis to model the stellar contribution in the Galac- 
tic extinction-corrected spectra in the rest frame (Cid Fer- 
nandes et al. 2005; Bian et al. 2007). We include 45 tem- 
plates from Bruzual & Chariot (2003; BC03) and one power- 
law component (representing the AGN continuum emission) . 
During the stellar population synthesis, we put twice the 
weight in the fit for strongest stellar absorption features, 
such as Call K, G-band, and Ca IIAA 8498, 8542, 8662 
triplet. For details, please see Bian et al. (2007) and the 
references therein. (2) The optical and ultraviolet Fe 11 tem- 
plate from the prototype NLS1 I ZW 1 is used to subtract 
the Fe n emission from the residual spectra after the above 
step. (3) Four Gaussians are used to model the H/3 profile, 
and two sets of two Gaussians are used to model the [O 
ill] AA4959, 5007 lines. We take the same linewidth for each 
component, and fix the flux ratio of [O in] A4959 to [O ill] 
A5007 to be 1:3. Two components of H/3 (from the NLR) 
are set to have the same linewidth of each component of [O 
ill] A5007 and their flux are constrained to be less than 1/2 
of each component of [O ill] A5007. Two broad components 
are used to model broad H/3 profile from the BLR contri- 
bution (broad component and intermediate component, BC 
and IC). We present our spectral fitting method in detail 
in Hu et al. (2008). In the above three steps, the best fit is 
reached by minimizing \ 2 , X 2 — ( Vi ~ y j aodo1 ) /N, where 
Oi is the error of the data set (xi,yi), N is the degrees of 
freedom. 

Objects without clear H/3 or [O ill] lines are eliminated. 
In order to obtain a reliable spectral fit, we carefully select 
objects for analysis. We select objects by the following cri- 
teria: (1) the signal-to-noise ratio (S/N) is larger than 15. 
The S/N is measured in the wavelength 4800-5040A, cov- 
ering the range for line-fitting. (2) \ 2 m above three steps 
(SSP, Fe 11 subtraction, and H/3 and [O ill] lines fitting) are 
less than 2.5. (3) the FWHM errors of BC and IC of H/3 from 
BLR and [O ill] AA4959, 5007 are less than 100%. (4) The 
height of one component of [O ill] is less than the height of 
one component of [O ill] by 30%. The first criterion leads to 
about 900 NLSls. The second and third criteria make sure 
we have a good fitting of H/3 and [O ill] . The fourth crite- 
rion make us reliable to obtain the core/narrow component 
of [O ill] line, and their line width is used to trace the cr*. 
Then we visually check these spectra one by one. In the end, 
we have a sample of 329 NLSls. 

The first moment of the line profile is 



A„ = 



/ XP(X)dX 
J P(\)d\ ' 



The second moment of the line profile is 
/ A 2 P(A)dA 



2 



JP(X)dX 



Xl 



(2) 



(3) 



We use the BC and IC to reconstruct the broad H/3 profile 
as -P(A). Then we use the method of Peterson et al. (2004) to 
measure the FWHM and ch/3 from the reconstructed broad 
H/3 profile (Equations (2-3)). The error of <th/3 is calculated 
from the error of FWHM for BC and IC in the line fit. The 
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Figure 1. Top: An example fit for SDSS J123831. 34+644356. 5. The black line is the original spectrum after Galactic-extinction correction 
in the rest frame. The red line is the contribution from the host galaxy. The green line is the Fe II emission. The residual is shown at the 
bottom. The left upper panel shows the region around Ca H+K AA 3969, 3934 and G-band. In the right upper panel we show the H/3 
and [O III] line fits, where the blue one with the highest peak is the IC, the blue one with the second highest peak is the BC, the two 
blue ones with the lower peaks are the components from NLRs. Bottom: An another example of line fit for SDSS J084716. 88+334858. 9. 
The black line is the original spectrum after Galactic-extinction correction, the starlight subtraction, and the Fe II subtraction in the 
rest frame. The bottom panel is the residual. The multiple Gaussian components are in blue and the sum of them is in red. The fitting 
window is in green. 



FWHM measurement is not sensitive to the broad wing of 
H/3, while the measurement of H/3 second moment is. We 
have to make sure the necessary of using BC and IC in 
the H/3 fitting (see also Dietrich et al. 2005). We use one 
broad component to model the H/3 profile from the BLR at 
the same time. When the x 2 in the fitting is decreased by 



20% with two-Gaussian respect to that with one-Gaussian, 
We have to make sure it is necessary to use BC and IC, 
and we use the result of two-Gaussian fitting, otherwise we 
use the result of one-Gaussian fitting. For H/3 profile with 
one-Gaussian fitting, the second moment is directly from 
our fitting, not from Equations (2-3). For all objects, the 
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measurement of FWHM is adopted from the two-Gaussian 
reconstructive profile. Hereafter, we call the total 329 NLSls 
as sample A, 209 NLSls with two-Gaussian fitting as sample 
B, and the other 120 NLSls with one-Gaussian fitting as 
sample C 

As the strongest forbidden line in the SDSS spectral 
wavelength coverage, we use the gas velocity dispersion of 
the narrow/core [O ill] component from NLRs to trace cr*, 

a foni] = V^obs [finst/(l + z)] 2 , where z is the redshift. 
For SDSS spectra, the mean value of instrumental resolution 
(Ti ns t is 60 km s" 1 for [O in] (e.g., Greene & Ho 2005a). 
We also reliably obtain er* (correction the resolutions of the 
SDSS spectra and BC03 templates) from our SSP synthesis 
for about 98 NLSls. The error of a, is given by different 
typical errors for different effective S/Ns at 4020 A (Bian et 
al. 2007). 



3 MASS FROM a H p AND FWHM 
3.1 Mass estimation 

With our measurement of o"h/3, we use the more recent 
Rblr — £5100 relation (which seems to hold for NLSls; Pe- 
terson et al. 2004) of Bentz et al. (2006) and / = 3.85 (Collin 
et al. 2006) to calculate the SMBHs masses in NLSls. We 
also calculate the mass from our FWHM for the H/3 profile 
from the BLR by Rblr — £5100 relation of Bentz et al. (2006) 
and / = 3 (Netzer 1990). In Figure 2, we compare these two 
masses. For sample B, the best linear fit with fixed slope of 
1 gives y=x-0.49, the mass from anp is on average larger 
by 0.49 dex than that from FWHM. In Figure 2a, we show 
the distribution of FWHM/cr, 1.33 ± 0.36, which deviates 
from 2.35 for a Gaussian profile. Our mass correction from 
ana respective to mass from FWHM is mainly due to the 
H/9 emission-line profile deviation from the Gaussian pro- 
file. For six NLSls in Peterson et al. (2004), we find that 
the mass based on a^p is on average larger by 0.46 dex with 
respect to that from FWHM, which is consistent with our 
calculation. 

By using the FWHM derived from the Lorentzian pro- 
file, Zhou et al. (2006) calculated SMBHs masses from the 
Rblr — L5100 relation of Kaspi et al. (2000, Rblr oc L°i 00 ) 
and / = 3. For comparison with the results of Zhou et al. 
(2006), we use the updated Rblr — £5100 relation of Bentz 
et al. (2006, Rblr oc I/5100 8 ) and / = 3 to calculate the 
mass. For the FWHM from Zhou et al. (2006), the updated 
Rblr — £5100 relation of Bentz et al. (2006) would lead the 
masses of Zhou et al. (2006) to be larger by 0.1-0.3 dex 
for NLSls with mass less than 1O 7 M0 (see Figure 2b). By 
the best linear fit through zero, we find FWHM in Zhou et 
al. (2006) derived from the Lorentzian profile is 0.84 ± 0.01 
of our FWHM for objects in sample B, which leads to the 
SMBH mass decreasing by 0.15 dex if we use the FWHM of 
Zhou et al. (2006). 

The usa of Rblr — £5100 relation of Kaspi et al. (2005) 
will decrease Mbh by 0.17 dex with respect to that of Kaspi 
et al. (2000). If we use / = 3.85 instead of 3, the mass will 
increase by 0.11 dex. The uncertainty of the mass calculation 
from the H/3 line is mainly from the systematic uncertainties, 
up to about 0.5 dex, which is due to the unknown kinematics 
and geometry in BLRs, and perhaps the effects of radiation 



Table 1. The SMBH masses from the a H p and FWHM. 





A 


B 


C 


log Mg H 

-^bol/-^Edd(°") 

AlogMg H 


7.53 ±0.50 
-0.81 ± 0.28 
0.03 ± 0.65 


7.69 ± 0.48 
-0.84 ±0.27 
0.12 ±0.68 


7.25 ±0.40 
-0.75 ±0.28 
-0.13 ±0.55 


log MfWHM 

iboi/i E dd(FWHM) 
AlogMfWHM 


7.17 ±0.44 
-0.44 ±0.37 
-0.34 ±0.65 


7.20 ± 0.44 
-0.35 ±0.37 
-0.38 ±0.70 


7.12 ±0.43 
-0.61 ±0.31 
-0.26 ±0.57 



f: A: for total 329 NLSls; B: for 209 NLSls in sample B; C: for 
120 NLSls in sample C. The mass deviation is based on Tremaine 
et al.'s relation when cr??T® is used for ct*. 



pressure (e.g., Krolik 2001; Peterson et al. 2004; Decarli et 
al. 2008; Marconi et al. 2008). 



3.2 The Mbh - of£ rc m] relation 

In Figure 3, we show the SMBH mass from (Th/s versus 
crfom]- I n Figure 3a, we show the mass derived from our 
FWHM versus c^q^. It is obvious that the mass from 
FWHM deviates from Tremaine et al. relation (solid line 
in Figure 3a), which is consistent with the result of Zhou et 
al. (2006). In Figure 3a, we show the SMBH mass from oh/3 
versus that from our FWHM. The SMBH masses based on 
(th/3 in 209 NLSls of sample B is larger by about 0.5 dex 
with respect to that from FWHM. 

For the solid line in Figure 3, the Mbh — cr* relation, 
M BH (cr*) = 10 8 - 13 [cr.7(200 kms" 1 )] 4 - 02 M (Tremaine et al. 
2002), we calculate the deviation of mass (see Table 1). In 
Figure 3b, we show the mass deviation from the Tremaine 
et al. relation versus the redshift. The dashed line is for no 
deviation, and the solid line is the relation found by Woo 
et al. (2008). For redshifts of NLSls larger than 0.4, their 
mean mass deviation is 0.28 ± 0.82 respect to Tremaine et 
al. relation, and for NLSls with redshift less than 0.4, the 
mean mass deviation is —0.04 ± 0.58 (Woo et al. 2008). We 
calculate the Eddington ratio, i.e., the ratio of the bolomet- 
ric luminosity (Lboi) to the Eddington luminosity (Lsdd), 
where L E dd = 1-26 x lO 38 (M B H/A/ )erg s _1 . The bolomet- 
ric luminosity is calculated from the monochromatic lumi- 
nosity at 5100A , Lboi = cbALa(5100A), where we adopt the 
correction factor cb of 9 (Kaspi et al. 2000; Richards et al. 
2006; Netzer & Trakhtenbrot 2007). We do not find larger 
mass deviations for objects with larger Eddington ratio. In 
Table 1, we give the mass, Eddington ratio, and the mass 
deviation from the Mbh — cr, relation for A, B, C samples. 
For sample B, mass from cths can be increased by ~ 0.5 dex 
with respect to that from FWHM. For sample C, mass from 
cr\ine can be increased by 0.13 dex with respect to that from 
FWHM, which is due to the values of / that we are using 
for sigma (/ = 3.85) and FWHM/2 (/ = 3) are not exactly 
consistent with a single Gaussian, whereas in order to get 
the same mass for a single Gaussian the ratio of f values is 
1.38. We also notice that, for NLSls with mass larger than 
10 8 Mq, they tend to fall above the Tremaine et al. relation. 
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Figure 2. The mass from the H/3 FWHM versus the mass from the H/3 second moment. The blue circles denote objects in sample B 
where the H/3 from BLRs can be fitted well by two-Gaussian with \ 2 decreased by at least 20% with respect to that by one Gaussian. 
The red circles denote objects in sample C with the H/3 from BLRs can be fitted well by one Gaussian in order to avoid the broad wing 
effect in measurement of the second moment. The solid line is 1:1, and the red line is the best linear fit with the fixed slope of 1 for sample 
B, y=x-0.49. Then FWHM and a^fi are measured from the reconstructed broad H/3 profile. For sample C, there is a scatter from 1:1. 
For sample C, mass from erii ne can be increased by 0.13 dex with respect to that from FWHM. (a) The distribution of FWHM/cr(H/3) 
for 209 NLSls in sample B. (b) SMBH mass from Zhou et al. (2006) versus that with updated the -Rblr — ^5100 relation of Bentz et al. 
(2006). 



4 DISCUSSION 

4.1 Large SMBH masses in NLSls? 

In Figure 3, we find that, for some NLSls, mass from cth/3 is 
very large, up to 10 s — 1O 9 M0, which is due to their obvious 
broad wings in their H/3 profiles. With the near-IR imag- 
ing data, Ryan et al. (2007) found that the average SMBH 
mass in their NLSls sample is 1O 7 ' 9 M0, where the mass is 
calculated from the mass and host galaxy luminosity rela- 
tion. The mass is typical for broad line Seyfert 1 galaxies. 
They find that it is larger by 1.5 dex respect to the mass 
calculated from the FWHM of Veron-Cetty et al. (2001). We 
note that these objects have broad wings in their H/3 pro- 
file, especially in their Ha profile (Figure 2 in Veron-Cetty 
et al. 2001). It is possible that the sample of Zhou et al. 
(2006) have some objects that should not be classified as 
NLSls. Considering the effects of the random velocity and 
the inclination in SMBH mass estimation of NLSls, Decarli 
et al. (2008) find that these effects would increase the mass 
for NLSls by 0.84 dex, which can account for the mass dif- 
ference between NLSls and BLSls. The definition of NLSls 
need to be reexamined if we think that NLSls harbor rapid 
growing small SMBH with high accretion rate. 



4.2 (j h/3 and FWHMs of BC, IC 

The flux ratio of H/3 to [O ill] A5007 from NLRs is often 
taken to be around 10% (e.g. Osterbrock & Pogge 1985; 
McGill et al. 2008). Rodriguez- Ardila et al. (2000) suggested 
this ratio is about 20%-100%, this high value due to that 
they just used two-Gaussian to model the total H/3 profile 
and assumed the narrow H/3 component is from NLRs. The 
flux ratio of H/3 to [O ill] A5007 from NLRs depend on the 
physics of the low-density gas found in NLRs photoionized 
by AGN, which is beyond the scope of this paper (e.g. Groves 
et al. 2006; Kewley et al. 2006). During our fitting procedure 
of the H/3 and [O ill] lines, we add a conservative constraint 
to the H/3 contribution from NLRs, less than a half of [O ill] 
A5007 flux. 

For 209 NLSls in sample B, the second moment <th/3 is 
calculated from the reconstructed H/3 profiles from BC and 
IC. We did a comparison between cth^ and the FWHM of 
IC/BC. In Figure 4, we show <t H /3 versus FWHM of IC and 
BC. The dashed line is the relation of FWHM = 2.35 x a 
when a Gaussian profile is measured. We find that there is a 
much stronger correlation between the an/3 versus FWHM 
of BC with respect to that for IC (see Figure 4). If using 
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Figure 3. The SMBH mass versus the core/narrow [O III] gas velocity dispersion. The mass is derived from the H/3 second moment. The 
sybols are same as Figure 2. The solid line is the Mbh — cr* relation of Tremaine et al. (2002). The green dashed line is the A/bh — cr* 
relation of Ferrarese & Ford (2005). (a) SMBH mass from H/3 FWHM versus the [O III] gas velocity dispersion, (b) Mass deviation versus 
the redshift. The solid line is the relation of Woo et al. (2008). 



the <7h/3 to calculate the SMBHs masses, the result depends 
much more on the BC FWHM, and less on the narrower IC 
FWHM. We think that if the H/3 profile from BLRs can be 
fitted well by one Gaussian, for mass, there is no difference 
for the usage of H/3 FWHM and the ct H /3. 



The existence of BC and IC in the H/3 profile from 
BLRs have been investigated by many people (Baldwin 
et al. 1998; Wills et al., 1993; Brotherton et al. 1994; 
Rodriguez-Ardila et al. 2000; Leighly 2004; Dietrich et al. 
2005), and it is suggested that BC and IC are emitted 
from two distinct emission region. For 12 NLSls, Dietrich 
et al. (2005) found that FWHMbc = 3275 ± 800km s" 1 and 
FWHMic = 1200±300km s -1 . For our 209 NLSls in sample 
B, we also found that FWHMbc = 4098 ± 1751km s" 1 and 
FWHMic = 1385 ± 492km s" 1 . Therefore, it is clear that 
there exists a BC in NLSls that are typical for BLSls, and 
the IC displays typical H/3 FWHM in NLSls. The equiva- 
lent width (EW) of BC for our 209 NLSls in sample B is 
19.4 ± 8.6A and 16.1 ± 8.0A for IC EW. There is no corre- 
lation between them. 



4.3 A/bh — <t * relation 

With FWHM from the H/3 Lorentzian profile, Zhou et al. 
(2006) use the -Rblr — £5100 relation of Kaspi et al. (2000) 
and / = 3 to calculate the mass. The updated -Rblr — £5100 
relation of Bentz et al. (2006) would lead to the mass of Zhou 
et al. (2006) larger by 0.1-0.3 dex for NLSls with mass less 
than 10 7 M Q (see Figure 2b), which would place NLS1S close 
to the Tremaine et al. Mbh — o", relation when cr* is adopted 
from ct[nii] (see Figure 29 in Zhou et al. 2006). For a sample 
of 58 NLSls selected from 11th edition of the "Catalogue of 
Quasars and AGNs " (Veron-Cetty & Veron 2003) and SDSS 
DR3, by the -Rblr — £5100 relation of Kaspi et al. (2005), 
Komossa & Xu (2007) suggested that NLSls do follow the 
Mbh — c* relation found in inactive galaxies after excluding 
"blue outliers" (see Bian & Zhao 2005). 

Using the second moment of broad H/3 profile, we find 
that the SMBH mass would be larger by ~ 0.5 dex with 
respect to that from H/3 FWHM in sample B with necessary 
two-Gaussian fitting. The larger masses make NLSls follow 
the Tremaine et al. relation when the cr^nii is used as a 
surrogate for the bulge velocity dispersion. In Figure 3, we 
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Figure 4. The FWHM of BC (top) and IC (bottom) versus the 
H/3 second moment for sample B of 209 NLSls. The dashed line is 
the relation between FWHM and the second moment, FWHM = 
o-jj/3 x 2.35, for a Gaussian line profile. 
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Figure 5. The Mbh — o", relation for 37 AGNs with reliable mea- 
surements of cr» . The blue and red circles have the same meaning 
as in Figure 2. The solid line is the Mbh — °* relation of Tremaine 
et al. (2002). The Green dashed line is the A/bh — °"» relation of 
Ferrarese & Ford (2005). The line in the left corner denotes the 
typical error in mass calculation. 



find that some AGN lie far below the A'Ibh — o~, relation and 
we tried to determine if they are "blue outliers". No "blue 
outliers" are found or it is impossible to measured the [O 
ill] blueshift due to the low S/N in [S n] or [O n] . We also 
use the [N n] FWHM (Zhou et al. 2006) as the a,, and find 
as similar result to that shown in Figure 3, but with more 
scatter. The [N n] FWHM measurement depends on the 
Ha profile fitting. And the [N n] FWHM is consistent with 
FWHM of [O in] core component, although the correlation 
is very weak (Komossa & Xu 2007). 

In Figure 5, we plot Mbh from au/3 versus <r« for 37 
NLSls. For the cr* from SSP synthesis, its uncertainty based 
on effective S/N at 4020A is typically about: 24 km s _1 at 
S/N=5; 12 km s" 1 at S/N=10; 8 km s -1 at S/N=15, where 
the effective S/N is the S/N (measured between 4010 and 
4060 A) multiplied by the stellar fraction (Cid Fernandes et 
al. 2005; Bian et al. 2007). The mass is calculated for AGNs 
that satisfy the first two criteria in section 2. The blue cir- 
cles denote AGN with necessary two-Gaussian fits and the 
red circles denote AGN with one-Gaussian fits, which is the 
same as that in Figure 2. For the total 37 NLSls in Figure 
5, the mass is between 10 6 — 10 8 Mq, and the distribution 
of the mass deviation is -0.24 ± 0.46. For 14 AGN with 
two-Gaussian fits, the distribution of the mass deviation is 
-0.14 ± 0.50, and for 23 AGN with one-Gaussian, the dis- 
tribution of the mass deviation is —0.29 ± 0.43. If these 23 
NLSls follow Tremaine et al. relation, we need / = 7.7. 
The ofoni] is slightly larger than a, for these 37 NLSls (see 
Botte et al. 2005). For about 3000 Seyfert 2 galaxies, Zhou 
et al. (2006) also found that <7[ NII ] — a, x (2.62/2.35) and 
overestimate the ct* (also see Onken et al. 2004). Therefore, 
using the gas velocity dispersion to trace ct» will place NLSls 
close to the Mbh — o"* relation and make the mass deviation 
smaller by 0.19 dex. 

For the total 37 NLSls, there exists marginal evidence 
that they deviate from the Mbh — c* relation found in in- 
active galaxies. For 23 NLSls with mass lower than 10 7 Mq, 
the deviation becomes much larger, up to — 0.46T0.40 (also 
see Figure 3), which is consistent with the result of Botte et 
al. (2005, their figure 3), although it is not the case for the 
sample of Greene & Ho (2006) (also see Barth et al. 2005). 
If these 23 NLSls follow the Tremaine et al. relation, we 
need /= 11.1. It is possible that the sample of Zhou et al. 
(2006) have some objects that can't be classified by NLSls. 
It is possible that there exists true NLSls with rapid grow- 
ing small SMBH. The reliable measurements of mass and 
the cr* are important for this kind of work. 



5 CONCLUSIONS 

The second moment of H/3 line from BLRs is calculated to 
derive the SMBHs masses for a sample of 329 NLSls selected 
from SDSS. The main conclusions can be summarized as 
follows: (1) For objects with necessary two-Gaussian fitting, 
the mean value of SMBH masses from the H/3 second mo- 
ment is larger by about 0.5 dex with respect to that from the 
H/3 FWHM. (2) Using the narrow/core [O ill] gas velocity 
dispersion as a surrogate for the stellar velocity dispersion, 
we find the new masses based on the H/3 broad emission line 
second moment bring them to the Tremaine et al. relation; 
(3) The H/3 second moment is more strongly correlated with 
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the BC FWHM rather than the IC FWHM. (4) Using the 
<r* measured from SSP synthesis, we find that, for NLSfs 
with masses lower than 10 7 Mq, they are marginally below 
the Tremaine et al. relation considering the larger scatter in 
mass calculation. If these 23 NLSls follow Tremaine's rela- 
tion, we need / = II. 1. 
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